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I The time-dependent variational approach to the pure Yang-Mills gauge theory, especially 

a color sit(3) gauge theory, is formulated in the functional Schrodinger picture with a Gaus- 
sian wave functional approximation. The equations of motion for the quantum gauge fields 
are formulated in the Liouville-von Neumann form. This variational approach is applied in 
order to derive the transport coefflcients, such as the shear viscosity, for the pure gluonic 
^ . matter by using the linear response theory. As a result, the contribution to the shear viscos- 

' ity from the quantum giuons is zero up to the lowest order of the coupling g in the quantum 

(~| ' gluonic matter. 

§1. Introduction 

> \ 

CN , One of recent interests for the quark and gluon physics, governed by the quan- 

' turn chromo dynamics (QCD), is to investigate properties of the quark-gluon plasma 

■ (QGP) and/or the quark-gluon matter. In the recent progress of the Relativistic 

I Heavy Ion Collider (RHIC) experiments, it is said that the QGP may be not free gas 

' but the strongly interacting quark-gluon matter.^) The matter composed of quarks 

Q!^ ■ and giuons seems to reveal the properties of the liquid, not gas, like the perfect liquid. 

This conjecture is derived by comaparing the obtained experimental data with the 
phenomenogical analysis by using the hydrodynamical simulation with rather small 
^ . shear viscosity, which leads to the near perfect liquid. The small shear viscosity is 

^ I also near the lower bound which is conjectured in the AdS/CFT correspondence.^^ 

Many works to understand properties of the gluonic matter were performed 
recently,^^"^^ while the transport coefficients, especially the shear viscosity for pure 
gluonic matter, were evaluated up to the lowest order of the QCD coupling constant 
g in the early study.^^ Namely, the shear viscosity rye for the gluonic matter at 
temperature T can be expressed as^^'^^ 

'nc = df— — , (1-1) 

r iog(i/r) 

up to the lowest order of g. Here, df is numerically determined constant. For the 
quark matter, the shear viscosity is also evaluated, for example, in the Nambu-Jona- 
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Lasinio (NJL) modeF^ by using the linear response theory,'^)' in which a rather 
smah shear viscosity is derived. ^''^ The smah shear viscosity leads to the short mean 
free path in general. Thus, it may be shown that the constituents, namely quarks 
and gluons, of the matter under consideration, are strongly correlated. 

However, as was shown by Ref.ll), the anomalous contribution to the shear 
viscosity in the turbulent plasma fields gives the small shear viscosity even in the 
weak coupling QCD in which the quarks and gluons are weakly correlated. Namely, 
the shear viscosity rj can be expressed by including the anomalous contribution 7]^ 
asii) 

ri~' = Vc'+VA ■ (1-2) 

If the anomalous viscosity tja is small, the total shear viscosity rj becomes to small 
even if the usual shear viscosity rjc is large for the small coupling constant g in 
Eq. (|l-ip . Thus, the small shear viscosity does not always lead to the strong coupling 
QCD, namely strongly correlated quark-gluon matter. 

In this paper, thus, we consider the pure gluonic matter as the weak coupling 
system, we are concentrated our interest to calculating the shear viscosity for pure 
quantum gluonic matter without quarks by using the linear response theory. In this 
paper, the quantum gluon means the quantum fluctuation part around the mean 
field which leads to the Eq. (|l-ip . Thus, the shear viscosity of the quantum gluonic 
matter gives the contribution of the next and higher order of g in comparison with 
Eq. (|l-l|) . One of purposes in this paper is to investigate the shear viscosity under 
small QCD coupling g for the quantum gluonic matter. It is important to investigate 
the contribution to the shear viscosity from the quantum gluonic fields. The reason 
is as follows: If there is a contribution to the shear viscosity of the order of g^, the 
finite value of the shear viscosity is remained even if the coupling g is small. 

To deal with the quantum gluons and to investigate the dynamics of the quantum 
gluons, the time-dependent variational method with the Gaussian functional as a trial 
wave functional in the functional Schrodinger picture may gives a useful tool. The 
reason why is that the mean fields and the quantum fluctuations around them can be 
treated on an equal footing and the higher order contributions for g are automatically 
included because certain kinds of the Feynman diagrams are taken into account in 
this variational approach. In this variational approach, the equations of motion 
for the mean fields and the fiuctuation modes around them are derived in a self- 
consistent manner. Especially, the equations of motion for the quantum gluon fields 
are formulated in a form of the Liouville-von Neumann equations. 

Another merit to use the time-dependent variational method for the pure gluonic 
matter is that the expectation values for various field operators and their products 
can easily be calculated because the state or the wave functional is prepared in the 
process of the variational calculations. When the transport coefficients such as the 
shear viscosity are calculated by using the linear response theory, the expectation 
values or thermal averages for the various operators such as the energy-momentum 
tensor operator are necessary. Thus, the variational approach may be suitable and 
give a practical method to calculate the transport coefficients. 

This paper is organized as follows. In the next section, the time-dependent 
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variational approach to the pure Yang-Mihs theory, especially the su{3) gauge theory 
as the QCD, is formulated in the Hamiltonian formalism. In §3, the time-dependent 
variational equations for the quantum gluon fields arc reformulated in a form of 
the Liouville-von Neumann equation for the reduced density matrix of the quantum 
fluctuation fields at zero and the finite temperatures. In §4, the shear viscosity is 
evaluated in our framework for pure gluonic matter by using the linear response 
theory from the viewpoint of weakly coupled QCD or weakly correlated pure gluonic 
matter. The last section is devoted to a summary and concluding remarks. 

§2. Time-dependent variational approach to QCD 

In this section, we give a variational method for the pure Yang-Mills gauge the- 
ory with color su{3) symmetry in the functional Schrodinger picture with a Gaussian 
approximation, which is developed in Ref.l2), in a slightly different manner. The 
trial state is constructed, paying an attention to the canonicity condition^^^' in 
our time-dependent variational approach. As a result, the equations of motion for 
variational functions are obtained as canonical equations of motion in classical me- 
chanics. 

2.1. Hamiltonian formalism of pure gauge theory 

In this subsection, we summarize the Hamiltonian formalism of the pure Yang- 
Mills gauge theory for the sake of the definiteness of notations. 

Let us start with the following Lagrangian density for the pure gauge theory 
with the color su{N) symmetry: 

F^, = d^Al - d^Al + gfabcA'^Al , (2-1) 

where A"^ represents the gauge field and the Greek indices such as n, u etc. and the 
Roman indices such as a, 6 etc. mean the Lorcntz and the color indices, respectively. 
The repeated indices are summed up. Later, we use another Roman indices such as 
i, j, ■ ■ ■ , which means the space components of the Lorentz indices, that is 1, 2 and 
3. Here, g represents the coupling constant and fabc is the structure constant for the 
color su{N): 

[Ta,n]= ifabcTc , (2-2) 

where {Ta} is the su{N) generators. In the adjoint representation, the su{N) gen- 
erator can be expressed as {Ta)l = —ifabc 

The conjugate momentum, tt"'*, for the field A^ is defined as 

dC 

tt"!" = = F"!" . (2-3) 
dA- 

Here, A = dA/dt. We introduce the vector notation such as = {A\,J^,J^). 
Then, the conjugate momentum with a space component can be expressed as 

_a / TTia TTia jpa \ 
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= -A'' -VA^o+9fabcA'A^o 

= . (2-4) 

Here, we define the color electric field E"". As is similar to the color electric fields, 
we define the color magnetic field as 

= V X A'^ - ^gfabcA" X . (2-5) 

Thus, we define the Hamiltonian density TCq as 
Ho = ^"^Al - C 

= i [{E'^f + [B-f] + 7r« • (V^g - gfabcA'Al) . (2-6) 

As is well known, the gauge theory leads to the constrained system. Namely, the 
conjugate momentum ttq is identical to zero, so it is necessary to impose a constraint 
condition and the consistency condition for the time evolution as 

7rg = F-^'O" = , ^0=0. (2-7) 

This fact leads to the Dirac theory of constrained system. In terms of the analytic 
mechanics, the constrained condition is written as 

< = { , j (fxH }p = £) . = , 

= V-E'' + ig- ifabcA' • E' , (2-8) 
where { , }p represents the Poisson bracket. Thus, the Hamiltonian is written as 

I d^xHo = I d^x^ [{E-f + (B'^f] + I d^xin'^ ■ V^g - gfabc^- ■ A'A^^) 

= I d'x^ [{E'^f + (B'^f] , (2-9) 

where we used the integrated by part in the second term and the constrained condi- 
tion ()2-7p . Thus, hereafter, we use the Hamiltonian density as 

no = l[{E''f + {B''f] ■ (2-10) 
For the later convenience, we introduce the following variable, Q: 

g = Q^T" 

y-E^ + ig- iUcA' ■ E'') r° 

= V -E + iglA', E'], (2-11) 

where A^ = A^^T'^ and so on, and Q"^ = D ■ tt". Thus, it is understood that Q is 
nothing but the infinitesimal generator of the gauge transformation. 
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2.2. Variational approach to pure gauge theory in quantum field theory 

In this subsection, we formulate the time-dependent variational method for the 
pure Yang-Mills gauge theory by using the functional Schrodinger picture^^^' 
within the Gaussian approximation. We formulate our variational method in the 
canonical form by the help of the canonical variable or canonicity conditions. ^^^'^^^ 

The time-dependent variational principle is formulated as 



<5 j dt{^\i-^ - J d^xH\^) = , (2-12) 

where 7i means the Hamiltonian density under consideration. In the functional 
Schrodinger picture, the commutation relation [ A1{x) , Ej{y) ] = i6ij5ab^^{x — y) 
leads to 

It is restricted ourselves that the trial state |^) or the trial wave functional 
<1>{A"') = (A"|(p) has the following Gaussian form as 

^{A'')=M-^eicp{i(E\A-A))e^^{-{A-A\^-iIJ\A-A)^ . (2-14) 

Here, we used abbreviated notations such as 

(E\A) = j d^xE''{x,t) ■ A^x) , 

^^'ib'^^ = // d'xd^y^ti^)lG-'^^{x,y,t)A';{y) . (2-15) 

Here, A'^{x,t), £'"(a;,t), Gf-{x,y,t) and Sfj{x,y,t) are the variational functions 
which are determined by the time-dependent variational principle. The reason 
why the form (j2-14p is adopted is that the canonicity conditions for {A^ , E^ ) and 
{Gf-,Z!fj) are automatically satisfied: 

{^\i±,\^) =£1 , (^|,^|<J)=0, 
6Ai bEi 

(<^l^^l^) = 0> = -Gl^ ■ (2-16) 

Thus, our time-dependent variational method is formulated as a canonical form. 

In the functional Schrodinger picture, the expectation values are easily calculated 
such as follows: 

{^A'iixm =X{x,t) , 

mEt{xm = Et{x,t) , 

{'!>\A'i{x)A]{ym = t)A]{y, t) + G'i^{x, y, t) , 

{^Et{x)E]{ym=W,{x,t)^{y,t) + \G-^'i^{x,y,t)+A{SGE)'^^ , 

{^A''{x)-E\x)\^) = A''{x,t)-E\x,t)+2{Gi:)'i^{x,x,t) . (2-17) 
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Thus, it is understood that represent the classical fields of gauge fields and 
the diagonal component of Gfj, that is, Gff , where indices i and a are no sum, 
is a quantum fluctuations around the classical field A". Thus, in this functional 
Schrodinger picture, the two-point function Gfj{x,y,t) plays a role of the gauge- 
particle propagator. 

It should be noted here that the trial state ()2-14p does not have the gauge 
symmetry, that is g\^) 7^ 0. Thus, we impose the gauge invariance by introducing 
the Lagrange multiplier. From (j2-8p . the constraint D -tv"" = is recast into another 
form Q°- = from (12-lip . where Q"" is the generator of the gauge transformation. 
Thus, we introduce the effective Hamiltonian density Ti by considering the gauge 
invariance in the space of the trial states as^'^^ 

n = Tio - uj" {x)g'' {x) , (2-18) 

where a;'* (a;) represents a Lagrange multiplier, which insure the constraint Q"'{x) = 0. 
Thus, we use the above Hamiltonian density in order to determine the time depen- 
dences of the variational functions A,- {x,t), {x,t), G'^j{x,y,t) and Z!fj{x,y,t). 

The expectation value of the Hamiltonian can be expressed as following simple 
form: 

{H) = (^1 1 d'x [Ho - uj''{x)g''{x)] \<P) 
= {Ho)- J d''xuj%x){g%x)) , 

(Fo) = j d''x{^-B\x)-B\x) + ]^E\x)-E\x) + ^Tr(^IG-V) 

1 ■ 9 
+2Ti{x\I:GS\x) + -Ti{x\KG\x) + — {TilS'T" {x\G\x)]) 

2 8 

+ ^Tr [S'T''{x\G\x)S'T''{x\G\x)] ^ , 

{g^{x)) = V -E^x) - igT,{x\T^[U , G\\x) + ig-ifabcA''-E\ (2-19) 
where we define 

Bi = eijkdjA,^ - -gf'^^'^eijkAjAf, , 
K = {-iS ■ Df - gS ■ B , 

D = V -igA , A = A^T^ , B = Wt"" . (2-20) 

Here, we can use the abbreviated notation such as {x\G\y) = Gfj{x,y,t). In the 
above representation, S represent the spin 1 matrices whose spatial component with 
i is S\ 

2.3. Variational equations and their solutions in the time-independent case 

The equations of motion for the variational functions are derived from the time- 
dependent variational principle in Eq. (|2T2p with the Hamiltonian density (j2T8p . 
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The results are summarized in the form of canonical equations of motion as 

-A{x,t)= , E\x,t) = - ^i^\ , (2.21a) 

6E {x,t) 6 A (x,t) 



'^r'^'^ix V t)- ^^^^ jj^hx V t)- ^^^^ 

6Z:fnx,y,t) 6Gf°{x,y,t) 



where {H) is given in Eq. (l2-19[) . 

In the time- independent case, the above equations of motion in Eq. (|2-21ap allow 
the following solutions within the lowest order of g as 

A" (a;) = , E'^ix) = -Vw"(a;) . (2-22) 

In the above solutions, the mean field or the classical field A'^(x) is identical to zero, 
namely, under this situation, only the quantum field is dealt with in this matter 
system. Further, in the lowest order of g, the time-independent equations ()2-21bp 
present the following solutions under A = : 

Gt^ik) = 5^' - ^) Gfc , Gk = ^ , (2-23a) 



d^k d^k' 
(27r)3 (27r)3' 

{k'a\U,,\bk) = \ [b,, - [b,^ - ^) gu^\q)b\]^ - k - q) 



uj^iq) = j d^xuj''{x)e-'^ '' . (2-23c) 

Thus, Gfj{k) has only the transverse component. This feature is plausible for the 
gauge-particle propagation. 

2.4. Thouless-Valatin correction 

In the Hamiltonian density (l2T8p . the constrained term uj"'Q"' is introduced. This 
treatment is resemble to that of the nuclear rotation. In the nuclear many-body 
theory, the collective rotational motion of the axially symmetric deformed nuclei for 
z-axis is described in the same way used in this section. If the nuclear rotation occurs 
in the perpendicular to the x axis, the state \^{t)) = e~^^^\^) is replaced as 

\^^{t)) = e~''^"-e~'^^^\'F) , (2-24) 

where the angular momentum operator and the angular velocity oj are introduced. 
Then, the Schrodinger equation, idt\\I^i^{t)) = H\<P^{t)) is recast into 

iH-Lol)\^) = E^\W) , (2-25) 
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where H is the original nuclear Hamiltonian. It is known that we have to get rid of 
the effect of nuclear rotation from the total energy as 

E = mH\^) - ^ ' = {^H\^) - AEty , (2-26) 

where / is the moment of inertia and is defined as^^-* 

/ = hm . (2.27) 

oj^O UJ 

This energy correction, AE'yy, in Eq. (j2-26p is well known as the Thouless-Valatin 
correction. 

Thus, in the approach to the pure Yang-Mills theory, it is first pointed out that 
the same correction term is necessary in Ref.l2). For the pure Yang-Mills theory in 
the treatment of the variational method, the Thouless-Valatin correction term can 
be expressed as 

AETy = j j d^xd^y{^\g^{x)g\y)\^){ax\^\by) , (2-28) 

where the moment of inertia for the gauge rotation, I"'^{x, y) = {ax\I\hy), is defined 
as 

T'^Hxy)- hm (2-29) 

It is first shown that, in Ref.l2), the above Thouless-Valatin correction term and the 
contribution of the moment of inertia for the gauge rotation play essential roles in 
order to reproduce the one-loop running coupling constant. The validity of our time- 
dependent variational approach owes the fact that the one-loop running coupling 
constant is exactly reproduced in the lowest order approximation of g under A = 
developed in Ref.l2). 

§3. Time-dependent variational equations for quantum gauge fields 

In this section, we present the equations of motion for the quantum fiuctua- 
tions around the classical field configurations A and E , namely, Gfj and E^j for 
quantum gauge fields, in a slightly different forms from Eq. ()2-21bp . We can formu- 
late the equations of motion for quantum gauge fields as the Liouville-von Neumann 
equation. 

3.1. Liouville-von Neumann equation for quantum gauge fields 

First, the reduced density matrix M is introduced as is similar to the Hartree- 
Bogoliubov theory for many-body physics in the boson systems. We define the 
reduced density matrix^'')' for the quantum gauge fields as 

j^ab. „ ,^ ^ ( -i{Af{x,t)El{y,t)) - i {Af{x,t)A>'^{y,t)) \ 

'^■^'^'^ I {En^,t)E^iy,t)) i{E-{x,t)A]{y,t))-\ 
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-2i{GE)'}^{x,y,t) Gf^{x,y,t) 
l{G'')t^{x, y, t) + A{UGS)1^{x, y, t) 2i{UG)f^{x, y, t) 

A1{x,t) = A'iix) - {A-ix)) = A-{x)-^{x,t) , 

Et{x,t) = Etix) - (Etix)) = Etix) - Et{x,t) , (3-1) 

where the symbol (• • • ) represents the expectation values for the state |<^(t))in (j2-14p 
as is shown in Eq. ()2-17p . In the later, the expectation values are replaced into the 
thermal averages. Here, Af means the quantum fluctuations around the classical 
configuration A^ . Thus, this reduced density matrix M can be regarded as the one 
consist by the quantum gauge fields. 

By the help of the Heisenberg equations of motion in the Heisenberg picture, the 
time evolution of the reduced density matrix is easily derived. As the result, we can 
obtain the following Liouville-von Neumann type equation of motion for the reduced 
density matrix composed of the quantum gauge fields as 

iMt^{x,y,t) = [n, M ]'}!;{x,y,t) , (3.2a) 



^2 



■HtHx, y,t)= "XL L ^aiL^ ) ^'(^ " ?/) , (3-2b) 



r^l'{x,t) = Kf^ + ^S^T'^x\G\x)S^T'^)f^ + |- {S^T'')%Tx \S,[T' l,x\G\x)\ 
\f^{x) = -igio'^ix)r% = gicu'^ixmt^ , (3-2c) 

where Kfj, and T'^ have been defined in Eq. (j2-20p . Here, Ti is the Hamiltonian 
matrix which governs the time evolution of the reduced density matrix. 

At this stage, it is important to indicate that the square of the reduced density 
matrix Ai is easily obtained from the second line of Eq. (j3-ip as 

M^^[l ?) , (3.3) 

Thus, the eigenvalues of the reduced density matrix are ±1/2 as is similar to the 
case of the linear sigma model. ^'^^'^^^ 

The eigenvector for the eigenvalue 1/2 can be expressed as 

{x\l/2nai) = (''rlt*}) , (3-4) 



Vnfix,t) 

where n represents a certain quantum number. Then the following eigenvalue equa- 
tion should be satisfied: 

From the above eigenvalue equation, we can derive the following equation: 

d'vM'^'lx V t) ( \---( ^ (3-6) 

d yM,, {X, y, t) ^ ^ J - 2 ^ J . (3 6) 
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Thus, we conclude that the v*) is the eigenvector for the reduced density 

matrix Ai with the eigenvalue —1/2, which is expressed as (a;| — l/2nai). Further, 
we can derive the following: 

Thus, we can introduce another vector as 

(x|l/2t«>=(<<^;5|) , (3.8) 

Then, we can express {l/2^nai\x) = (f *f (a;, t), u*f (cc, t)). 

3.2. Spectral decomposition of reduced density matrix 

In the previous subsection, it is learned that the reduced density matrix M has 
the eigenvalues ±1/2. In this subsection, the eigenvalue equation is summarized 
taking into account the extension to the finite temperature systems. Further, by 
using the eigenstates for 7W, the reduced density matrix is expressed in the form of 
the spectral decomposition. 

Considering the extension to finite temperature systems, the eigenvalue equa- 
tions are described with the abstract representation as 

Mt^\afnA3) = <yfn\ofn,a,i) , Mf\aflb,j) = af„\afla,i) , (3-9) 

where / = 1/2 and cr = it at zero temperature developed in the previous subsection. 
From the second equation in (j3-9p . we obtain 



{afl,b,j\M'^'i = aU{afla,i\ . (3-10) 
By using Eqs. (|3-9p and (j3-10p . we can easily derive the following orthogonal relations 



as 



Y{(^'fl',a,i\(jfn,a,i) = , (3T1) 

a,i 

where the normalization condition of \afn,a,i) is taken into account. By using 
the completeness relation J d"^x\x){x\ = 1 and the expression such as (|3-4p with 
c7 = cj' = loro" = l and a' = —1, the above condition (j3Tip is rewritten as^^^^ 

Yl j K't{x)Unt{x) + ul,1{x)Vnt{x)) = b^n' , 

Y I iu:,t{x)v:tix) - v:,t{x)uinx)) = o , 

Y J d^'x {un'1{x)vnt{x) - (a;)n„^(a;)) = . (3-12) 
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Thus, the reduced density matrix M. itself can be expressed in terms of the eigen- 
states of A4 as follows: 



Mt^{x,y,t)= ^ U 

n(o->0) 



+^ \ . .1 

(3-13) 



We can easily verify that the above M satisfy the eigenvalue equations (j3-5p and 
(j3-6p with fn instead of 1/2. 

Next, let us determine the eigenvalue in the system at finite temperature. In 
the time- independent case, from Eq. (j3-2ap . it is seen that the reduced density matrix 
and the Hamiltonian matrix commute each other. Thus, there exist the simultaneous 
eigenstates for M. and TL whose eigenvalues are afn and a En, respectively. Namely, 

hi 



,* bi 

nj\ 
*b 

Thus, we can derive TT{7iM.) = Here, the Helmholtz free energy is defined 

as 

F = {n)-TS , 

{n) = Tv{nM) = Y,Em , 

m 

f^ma) f^ma^^T^ma) j (3T5) 

where T is temperature. Thus, the minimization condition is imposed: 

5F = TT{mM) - T5S = . (3-16) 

Here, we assume that afm depends on 72^0- linearly. Under this assumption, 
5M/5nma = v){v* ,u*) for cr = + and bM-jbumu = ^{u* , —v*){—v, u) for a = — 
are obtained respectively. Thus, we obtain 

5F 1 + nrnn 



j- -"lit 

I.e., Uma = g£;„/T _ ^ ' (^'l"^) 

where o" = it. Thus, we omit the suffix a in Uma- Of course, the eigenvalue / is 
reduced to 1/2 when T ^ 0. Thus, finally, we obtain 

fm — f^m + 7: ; 'T-m = 77^ ~ ■ (3T8) 

2 e^™/-' — 1 
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From the (l,2)-component of the reduced density matrix in Eq. (l3-ip . the two- 
point function G can be decomposed by each spectral function as 

Gt^{x,y,t) = Y,fn [unnx)<'j{y) + <t{x)u4{y)] . (3-19) 

n 

Hereafter, we take the quantum number n as momentum k. Then, we obtain 

Gt^{x,y,t)= I ^f^[ukt{x)ul'^{y) +uinx>kj(.y)] ■ (3-20) 

At zero temperature, we have already derived the expression of G in Eq. (j2-23p . In 
the above expression in (j3-2Up . = 1/2 and uu* is also expressed. By using the 
knowledge of the zero temperature case, the above expression in ()3-20p can be recast 
into 



G-(.,^,t) = 5-/ ^Gl (6., - ^) e-(--) , 

G^ = (2nfc + 1)^ . (3-21) 

From the eigenvalue equation (j3-14p . the energy eigenvalue for the quantum 
gauge fields is easily obtained. For simplicity, we neglect the gauge rotating term, 
namely, we put uj'^{x) = 0. Then from (j3T4p . we obtain 

rV2«^(fc) = 5'^^|fc|(<^.,-^) (3-22) 

up to the order of g. 

Finally, it should be noted that the gauge-particle has only transverse component 
as is realized in Eq. (j2-23p with factor {6ij — kikj/k^) . Since the reduced density 
matrix ^A includes G and E composed by the quantum gauge field, ^A also contains 
the factor {5ij — kikj/k^). Thus, we define the following projection operator P: 

Then, V certainly has a property of the projection operator, namely, "P^ = V. As is 
easily shown, the following relations are satisfied: 

G = VGV , S = PEV , 

M = VMV = VM = MV . (3-24) 



Since there always exists the projection factor to the transverse component, (5; 
kikj/k^), the inverse of the two-point function G should be regarded as the two point 
function which satisfies the following relation: 

G-^G = V. (3-25) 
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§4. Transport coefficients of quantum gluonic matter 

Hereafter, wc deal with the color su(3) pure gauge theory, namely the QCD 
without quarks. In this section, we present the expression of the transport coefficients 
in our variational approach by using the Kubo formula^)' for the quantum gluonic 
matter with the color su{3) symmetry. 

4.1. Kubo formula based on the variational approach 

By taking into account an external source field A and its conjugate force F{t), 
the Hamiltonian H is modified from Hq to 

H{t) = Ho- AF{t) . (4-1) 

If the external force F{t) is adopted as F{t) = Fe~^*, then, an observable 13 at time 
t,{B)t, can be expressed in the linear response theory as 

{B)t = (^)eq + XBA{u^)Fe-"^' , (4-2) 

where (• • • )eq means the thermal average with respect to the equilibrium state. Here, 
Xba{^^) is called the complex admittance and is defined as 

oo 

iuit—et 



Xba{co) = - hm - / dt{[A, B{t) ] )eqe 



oo 

hm / dtipBA{t)e'^'''' , (4-3) 



where 



^BA{t) = -t{[A, B{t) ] )eq 



/ d\{A{-ih\)B{t))^^ (4-4) 



is a quantum response function. By using the integrated by part and an property of 
the equilibrium state, namely, {AB)^^^ = —{AB)^^, the complex admittance is recast 
into 



1 r'^ 



XBAico) = — hm / dte-*-^*( [ B{t) , A(0) ] )eq 

1 r°° - 

-J^ i^oio ^ ' ^^^^ ^ ■ ^^'^^ 

Following the general theory,^-* the transport coefficients are obtained by adopt- 
ing both operators A and B being currents J{t) as 

XBAiu;,k = 0) = - hm Hdt [ d^re'^^'-^'i [ J{r,t) , J(0,0) ] )eq 

\im^ J^dt J d^re-''{ [ J{r,t) , J(0,0) ] )eq , (4-6) 
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where we return to the natural unit, h = 1. Here, current J{t) is defined as 

J{t) = J d^rJ{r,t)e-'^-'- (4-7) 

for an application to the gluonic matter in mind. 

4.2. Shear viscosity in the quantum gluonic matter 

In the gluonic matter, the dependence of the couphng constant g for the usual, 
not anomalous, shear viscosity is given as p-ip in the lowest order of gJ'^'^^ In this 
paper, since we deal with the quantum gluonic field described by the variables G 
and X! with A = 0, so we can evaluate the usual shear viscosity rjc with higher order 
of g compared with those developed in the previous papers. Hereafter, we denote 
r]c as r] simply because, as a final result in this paper, the small value of r]c can be 
derived from the viewpoint of the weak coupling QCD. Thus, it is expected that ry 
is also small from Eq. (jl-2p for the quantum gluonic matter. 

In order to calculate the shear viscosity, the energy-momentum tensor for the 
pure gluonic field is necessary. In the symmetric representation, the energy-momentum 
tensor Ty is obtained as 



Then, Tif can be expressed in terms of the color electric and color magnetic fields as 

Too(r) = ^(E'^ir) • £;"(r) + B-{r) ■ B\r)) , 
Toi(r) = -eijkE^{r)Bl{r) , 

T,j{r) = -E^{r)E^{r) - B^{r)B^{r) + hijiE'^ir) ■ E^{r) + S'^(r) • S'^(r)) . 

(4-9) 

The shear viscosity r/(a;) is obtained by taking the current J as T^y = T12 in 
Eq. (l¥6]l :^°) 



7]{u:) = - \n^{uj) - 

UJ '- 

771^(0;) = -i lim / dt / dVe^"*-^*( [ ri2(r,t) , ri2(0,0) ] )eq . (4-10) 

By taking a limit — > 0, the shear viscosity r/(0) has a simple form^''^ as 

d 



ry(0) = — ;:-im n^iuj) 
duj 



(4-11) 



Here, the thermal average (• • • )eq can be replaced to the expectation value given in 
Eq. dMTD for G and (l2^m]) for E with instead of Gk in our variational approach 
at finite temperature. Thus, it is necessary for calculating the shear viscosity to 
evaluate the thermal average {[Ti2{r,t) , Ti2(0,0)]) in our framework. 
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However, we need the operator at time t, namely, Ti2{r,t). In order to derive 
this operator at t, we only need to evaluate the operator E"" and B'^ at time t because 
the expression of Ti2{r, t) has same dependence with respect to E"' and B", that is. 



Ti2ir,t) = -Et{r,t)E^{r,t) - Bl{r,t)B^{r,t) 



(4-12) 



The operators E°'{r,t) and B°'{r,t) are obtained as a result of the time evolution 
governed by the Hamiltonian J cL^xTLq with (j2-10p : 



(r, t) = e^^o*£;f (r)e-*^«* , 5f (r, t) = e*^o*5f (r)e-^-^o* 
Ho = j d^rno = I d3ri[(£;'^(r))2 + (S'^(r))2] . 



(4-13) 
(4-14) 



Here, we can derive 
[ ^0 , Efir) ] 



-i{S-pUBt{r) + 0{g) , 



[ H, , Bt{r) ] = ie,,,dkE^{r) + -ge^,kr'^{E'^{r)AUr) + AUr)E^ir)) 



i{S-p\uEl{r)+0{g), 



(4-15) 



P 



dr ' 



where S is defined in (j2-20p . Thus, up to the lowest order of g in our quantum 
gluonic matter, we can derive the operators at time t as 



Etir, t) = Ct^(t)E'A^) + Dt^{t)B^ir) + 0{g) , 
Bt{r,t) = Ctf{t)B]{r) - D-^{t)E^{r) + 0{g) , 
C^j'it) = 6^'[cos{t{S ■ , Dt^it) = 6'^''[sm{t{S ■ p))] 



(4-16) 



Thus, we can derive the energy- momentum tensor operator at time t in Eq. ()4T2p . 

From (|4-10p . we need the ([Ti2(r, t), ri2(0, 0)]) to estimate the shear viscosity 
in quantum gluonic matter. After lengthy but straightforward calculation, we can 
derive the following form up to the lowest order of g: 

{[Tuir,t), ri2(r',0)]) 



e,u^^2\r,r',t) + e,2iS^,\ry,t) 



Nti,{t)-idl5[r-r') 



+e,nSf^{r',r,t) + e,2iSt^{r\r,t) 
e,iir,fir,r',t) + e,2iViir,r',t) 



eaiTiUr',r,t)-ei2ir^'{r',r,t) 



+0{g^) , (4-17) 
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where we define 



Mfj.jit) = cmc%{t) + DmD%{t) , 



Sf^iry, t) = e,^kd:^{GS)t^ (r, r', t) + e,^kdt{SG)t{r, r' , t) + 0{g^) 



cb , 



"icb , 



T^f{r,r',t) 



(5 • pMGE)tJr,r',t) - {SG)t{ry ,t){S ■ + 0{g 



ab I 



'^(G* ^)ij{'r',r',t) + eipm^jkndjjdl Gfl^^{r,r',t) 
-A{UGU)t^{r,r',t) + 0{g^) 

kG-')tHr, r', t) + {S- pUGtnir, r\ t){S ■ $)nj + 0{g 



(4-18) 



where we define f{r,r',t)'p = —idf{r,r',t)/dr. Here, we used the fact that U is 
of order of g as is seen in Eq. ()2-23p . 

4.3. Lowest order approximation for shear viscosity in quantum gluonic matter 

In the lowest order of g, M^2ij Eq. (14-18p is proportional to 6"'^. On the other 
hand, S^j is proportional to f°-'"^ui'^ because both the variables (Gi:)g and {IJG)^^ 
in Sfj are proportional to f°-^'^u)'^ which is derived from (l2-23bp with and ()3-2ip . 
Thus, in Eq. (j4-17p . the first term vanishes, namely. 



^ijki-mn = for any i, j,k, I, m and n . 



(4-19) 



Next, let us consider the second term in (j4-17p . At zero temperature, the equa- 
tion of motion for G can be derived from (j2-19p or ()2-21bp by 5{H)/6G = which 
leads to 



8 2 

in the lowest order of g. Thus, the solution of G is written as 

1 1 



G 



2Vk 2S-P 



(4-20) 



(4-21) 



where K = {S ■ p)^ with A = in (|2-20p . Using the above fact, from ()4-18p . the 
following is derived up to the lowest order of g: 



y^ab 



\g'^-{s-p)G{s-p) 



ab 







(4-22) 



As a result, the shear viscosity for the pure quantum gluonic matter is zero at zero 
temperature up to the order of g^ : 



r]{uj) = 



(4-23) 
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up to the order of g at zero temperature. 

At finite temperature, we have the solution for G and U in Eq. (j3-2ip and (|2-23bp 
with instead of Gk- Thus, we obtain T in the lowest order approximation as 



I,^[r,r,t)-d J 2\k\ ^d,, - 



' i^-)' sinhfl^U " ^ ^ 



i lim^ / dt I d'^r e'^^-'^ N'^!^^j{t){idi5\r)) 



T ^ 

From Eq. (j4-10p . the complex admittance for the shear viscosity r]{oj) is written as 
n^iio) = -i lim / dt / d'^r e'"*-^*([ Ti2{r,t) , ri2(0,0) ]) 

X [e,i,2;f (r, 0, t) + ej2iT,f{r, 0, t) - eaiT^^iO, r, t) - ei2iT^^iO, r, t)] 

= i lim^ / '^M e''^^-'^N^ijit)[6\r) 

xUjiiTtUr, 0, t) + e,2iTi%{r, 0, t) + eaiTi%{0, r, t) + e,2iri%{0, r, t)}] , 

(4-25) 

where the integration by part has been carried out from the second line to the third 
line and we define T} jj(r,r',t) as 

Y, ..(r r' t) - 6"^ I ik-(r-r') ki\k\ / 

''^^^ ' ' ^" y (2vr)3^ sinh(|^)r^ I^IV' ^ ^ 

Here, up to the lowest order of g in this quantum gluonic matter, the operator 
Ni2ij{t) is written in Eg. diTSl) with (Ii46|) . Thus, we can expand Nuijit) as 



12ii 



n=l 



= [6utiS ■ p)2j - 52jt{S ■ p)u + 0(p')] , (4-27) 

where f{xys are expansion matrices with numerical factors for this expansion and 
p = —id /dr. Thus, 

I d3riVi%(t) [6Hr)Ti:Ur,0,t)] = (4-28) 

for any i, j, k and m because the integrand is the total derivative with respect to r. 
Thus, we can obtain the relation 7I^(a;) = 0. Finally, from Eqs. (j4-25p and ()4-28p . 
we conclude that the shear viscosity rj{uj) in Eq. (|4-lUp is as follows 

r]{uj) = (4-29) 
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up to the lowest order of g in the quantum gluonic matter. Here, the terms of the 
next order of C'^j{t) and Dfj{t) in N^2-{t) include the couphng constant g without 
the spatial derivative. Thus, Eq. ()4-28p is not satisfied in the order of g. Therefore, 
the result (j4-29p is valid up to the order of g^ while (j4-23p is valid up to the order of 
g^ at zero temperature because the mechanism to vanish the value of shear viscosity 
is different. 

§5. Summary and concluding remarks 

In this paper, the time-dependent variational method for the pure Yang-Mills 
gauge theory is formulated in the functional Schrodinger picture with the Gaussian 
trial wave functional. In this variational method, the classical mean fields and the 
quantum fluctuations around them are treated self-consistently and both degrees of 
freedom are coupled each other. Further, the equations of motion for the quantum 
fluctuations around the mean fields were reformulated in a form of the Liouville- 
von Neumann equation for the reduced density matrix which was introduced in the 
Hartree-Bogoliubov approximation developed in the many-body problems for boson 
systems. 

This variational method developed in this paper was applied to the pure quan- 
tum gluonic matter system in order to evaluate the shear viscosity, which is one of 
the transport coefficients of the gluonic matter in the system with the color su{2>) 
symmetry, namely, the QCD without quarks. As a result, it was shown that there 
is no contribution of the quantum gluons to the shear viscosity in the pure gluonic 
matter up to the lowest order of the QCD coupling g at finite temperature. Namely, 
up to the order of g^, the contribution of the quantum gluons to the shear viscosity 
is nothing. At zero temperature, adding to the order of g^ , there is no contribution 
up to the order of g due to the equations of motion. Thus, for small g, namely, from 
the viewpoint of the weak coupling QCD,^^) the shear viscosity in quantum gluonic 
matter may be small because the quantum gluons contribute to the shear viscosity 
from the order of g or higher at finite temperature. 

Recently, Matsui and Matsuo give the transport equations for the Wigner distri- 
bution function and the anomalous distribution function^^-* which are coupled each 
other to determine the dynamics of the meson fields and the fluctuations around 
them in the linear sigma model, as is similar to our formalism. In order to compare 
the theoretical analysis with the experimental results, the information of the gluon 
distribution function may be necessary as was discussed in Ref.22) in the case of the 
linear sigma model. It is one of important further problems to investigate the gluon 
distribution function governed by the transport equation derived by this formalism 
such as the extended Boltzmann equations. 

It may be interesting to investigate the higher order contribution to the shear 
viscosity because the approximation used in this paper corresponds to the Hartree- 
Bogoliubov like approximation. The random phase approximation (RPA) is missing 
in the treatment in this paper. It may be necessary to extend our treatment to in- 
cluding the RPA like modes as was developed for the linear sigma model. Further, 
it may be also interesting to investigate the behavior of other transport coefficients 
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in this framework developed in this paper. They are future problems. 
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